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Weak lensing of CMB anisotropies and polarization for the power spectra and higher order statistics 
can be handled directly in harmonic-space without recourse to real-space correlation functions. For 
the power spectra, this approach not only simplifies the calculations but is also readily generalized 
from the usual flat-sky approximation to the exact all-sky form by replacing Fourier harmonics with 
spherical harmonics. Counterintuitively, due to the nonlinear nature of the effect, errors in the 
flat-sky approximation do not improve on smaller scales. They remain at the 10% level through 
the acoustic regime and are sufliciently large to merit adoption of the all-sky formalism. For the 
bispectra, a cosmic variance limited detection of the correlation with secondary anisotropies has an 
order of magnitude greater signal-to-noise for combinations involving magnetic parity polarization 
than those involving the temperature alone. Detection of these bispectra will however be severely 
noise and foreground limited even with the Planck satellite, leaving room for improvement with 
higher sensitivity experiments. We also provide a general study of the correspondence between flat 
and all sky potentials, deflection angles, convergence and shear for the power spectra and bispectra. 



I. INTRODUCTION 

As the cosmic microwave background (CMB) photons 
propagate from the last scattering surface through in- 
tervening large-scale structure, they are gravitationally 
lensed. Weak lensing effects on the the temperature and 
polarization distributions of the cosmic microwave back- 
ground is already a well-studied field. As in other as- 
pects of the field, early work treating the effects on the 
temperature correlation function ^ has largely been su- 
perceded by harmonic space power spectrum analyses in 
the post-COBE era In harmonic space, the physical 
processes of anisotropy formation are most directly man- 
ifest. However for weak lensing in the CMB, correlation 
function underpinnings have typically remained, forcing 
transformations between real and Fourier space to de- 
fine the effect in a small-angle (flat-sky) approximation. 
Exceptions include recent work on the non-Gaussianity 
of the lensed temperature field where a direct harmonic 
space approach has been taken 

In this paper, we provide a complete framework for the 
study of lensing effects in the temperature and polariza- 
tion fields directly in harmonic space. Not only does this 
greatly simplify the power spectrum calculations but it 
also establishes a clear link between weak lensing power 
spectrum observables in wide-field galaxy surveys and 
CMB observables for cross-correlation studies. Further- 
more, this approach is easily generalized to lensing on 
the full sky by replacing Fourier harmonics with spheri- 
cal harmonics. 

We show that counterintuitively, corrections from em- 
ploying an exact all-sky treatment are not confined to 
large angles. The second order nature of the effect brings 
in large scale power through mode coupling. Since the 
all-sky expressions are as simple to evaluate as their flat- 
sky approximations, which themselves are much simpler 
to evaluate than the correlation function analogues, they 



should be employed where full accuracy is required, e.g. 
for the analysis of precise measurements from CMB satel- 
lite missions. 

Beyond the power spectrum, lensing induces three 
point correlations in the CMB through its correlation 
with secondary anisotropies |^,^, even when the intrin- 
sic distribution at last scattering is Gaussian. Detec- 
tion of these effects in the temperature maps however 
are severely limited by cosmic variance. The primary 
anisotropies themselves act at as Gaussian noise for these 
purposes. In this case, the low level at which the CMB 
is polarized can be an asset not a liability. Three point 
correlations involving the polarization, where orientation 
plays a role, are most simply considered with their har- 
monic space analogue, the bispectrum. We introduce 
polarization and polarization-temperature bispectra and 
show that they can have signal-to-noise advantages over 
those involving the temperature alone. 

The outline of the paper is as follows. In §||, we treat 
the basic elements of the cosmological framework, CMB 
temperature and polarization, and weak lensing needed 
to understand these effects. Detailed derivations are pre- 
sented in a series of Appendices: ^ covers the all-sky 
weak lensing approach, |b| the evaluation of the all-sky 
formulae and ^ the correspondence between the flat and 
all sky approaches for scalar, vector and tensor fields on 
the sky. The lensing effects on the power spectrum are 
treated in the flat-sky approximation in §111 and in the 
exact all-sky approach in §IV. In §0 we study the ef- 
fects of lensing on the bispectra of the temperature and 
polarization distributions. We conclude in § [v| . 



II. FORMALISM 

In this section, we review and develop the formalism 
necessary for calculating lensing effects in the CMB. We 
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review the relevant properties of the adiabatic cold dark 
matter (CDM) model in §11 A. In §11 B, we discuss the 



power spectra and bispectra of the temperature fluctu- 
ations, polarization and temperature-polarization cross 



correlation. Finally in §IIC, we review the properties of 



weak lensing relevant for the CMB calculation. 



A. Cosmological Model 

We work in the context of the adiabatic CDM family of 
models, where structure forms through the gravitational 
instability of the CDM in a background Friedmann- 
Robertson- Walker metric. In units of the critical density 
3H^/8nG, where Hq = lOOh km s^^ Mpc^^ is the Hub- 
ble parameter today, the contribution of each component 
is denoted fij, z = c for the CDM, b for the baryons, A for 
the cosmological constant. It is convenient to define the 
auxiliary quantities = + ^6 and fix = 1 — X]j 
which represent the matter density and the contribution 
of spatial curvature to the expansion rate respectively. 
The expansion rate 



(1) 



then determines the comoving conformal distance to red- 
shift z. 
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FIG. 1. Temperature and temperature-polarization cross 
power spectra. Shown here are the power spectra of the un- 
lensed and lensed fields, their difference in the all-sky and 
flat-sky calculations and the error induced by using the flat 
sky expressions. The oscillatory nature of the difference indi- 
cates that lensing smooths the power spectrum. 



D{z) = 



H{z') 



-dz' , 



(2) 



in units of the Hubble distance today Hq ^ = 
2997.9ft,~^Mpc. The comoving angular diameter distance 



Da = n^'"^ sin\i{n'^^D 



1/2, 



(3) 
-0, 



plays an important role in lensing. Note that as VLk 

The adiabatic CDM model possesses a power spectrum 
of fluctuations in the gravitational potential $ 



k 



T\k), (4) 



where the the transfer function is normalized to T(0) = 1. 
We employ the CMBFast code ^ to determine T{k) at 
intermediate scales and extend it to small scales using 
analytic fits 0. 

The cosmological Poisson equation relates the power 
spectra of the potential and density perturbations 5 



Here 5h is the amplitude of present-day density fiuctua- 
tions at the Hubble scale; we adopt the COBE normal- 
ization for 5h |§- F{z)/(\ + z) is the growth rate of 
linear density perturbations 5{z) = F(z)(5(0)/(1 + z) [|j 



F{z) cx (1 z) 



H(z) 

Hn 



dz'{l + z') 



Hn 



H{z') 



(7) 



For the matter dominated regime where cx (1 -t- z)"^^^, 
F is independent of redshift. 

Although we maintain generality in all derivations, we 
illustrate our results with a ACDM model. The parame- 
ters for this model are Vic = 0.30, fif, = 0.05, I^a = 0.65, 
h = 0.65, Yp = 0.24, n = 1, and 6h = 4.2 x 10"^ This 
model has mass fiuctuations on the 8h Mpc~^ scale in 
accord with the abundance of galaxy clusters erg — 0.86. 
A reasonable value here is important since the lensing 
calculation is second order. 



B. CMB 
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nlil + zfAj, (5) 



and gives the relationship between their relative normal- 
ization 



AM = I 



4" 



K 



nlF{z)5l. 
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We decompose the CMB temperature perturbation on 
the sky 8(n) = Ar(n)/T into its multipole moments 



e(n) = ^e,„,y/"(n). 



(8) 



Im 



The polarization on the sky is represented by the trace- 
free symmetric Stokes matrix on the sky 
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FIG. 2. Polarization power spectra. The same as in Fig. |l| 
except for the E and B polarization. We have assumed that 
the unlensed B spectrum vanishes as appropriate for scalar 
perturbations. 



P(n) = +X{n) (m+ 
where 



_X(n) (m_ ® m_) , (9) 



-_X{n) = Q{h)±tU{n), 
m± = — (ee ip ie^,) . 



(10) 



The complex Stokes parameter ±X is a spin-2 object 
which can be decomposed in the spin-spherical harmonics 

m 



±X(n)-^±Xz„±2y,"(n). 



(11) 



We have assumed that the Stokes V parameter vanishes 
as appropriate for cosmological perturbations; for a full 
set add the term Veij to the polarization matrix, where 
eij is the Levi-Civita tensor. 

Due to the parity properties of the spin-spherical har- 
monics 



one introduces the parity eigenstates 

±Xlm — Elm ± iBim , 



(12) 



(13) 



such that Elm just like 8;,„ has parity (—1)' ("electric" 
parity) whereas Bim has parity (—1)'+^ ("magnetic" par- 
ity). Density (scalar) fluctuations in linear theory only 
stimulate the E component of polarization. 



The power spectra and cross correlation of these quan- 
tities is defined as 



XX' 



(14) 



where X and X' can take on the values ld,E,B. Note 
that the cross power spectra between B and 8 or i5 have 
odd total parity and thus vanish assuming anisotropy 
formation is a parity invariant process. 
The bispectrum is defined as 
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{XlmXiim' Xi,, m") = in 



, (15) 



and vanishes if the fluctuations are Gaussian. Even in the 
presence of non- Gaussian! ty due to non-linear but parity- 
conserving sources, bispectra involving an even number 
of magnetic parity terms (including zero) vanish for I + 
I' + I" =odd and those involving an odd number vanish 
for l + l' + I" =even. 

For a small section of the sky or high multipole mo- 
ments, it is sufficient to treat the sky as flat. In the 
flat-sky approximation, the Fourier moments of the tem- 
perature fluctuations are given as 



e(n) = 



{2nf 

and the polarization as 
dH 



e(l)e' 



,X(n) 



,X(l)e 



±2i(ipj-i/3)gil-n 



(16) 



(17) 



where ipi is azimuthal angle of 1. Again one separates the 
Stokes moments as 



±X{\) = E[\)±iB{\). 



(18) 



As in the all-sky case, the power spectra and cross 
correlations can be defined as with power spectra 



(X*(1)X'(1')) ={2irf5{\-\')Cf,) 



XX' 



(19) 



(X*(1)X'(1')X"(1")) = (27r)25(l - 1' - l")S(f.0)" 

The power spectra for the fiducial ACDM model are 
shown in Figs, p] and ||. 

In Appendix^, we establish the correspondence be- 
tween the all-sky and flat-sky spectra. For the power 
spectra and bispectra. 
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for sufficiently high Vs. 

For the bispectra, we have assumed that the triplet 
is composed of an even number of magnetic parity {B) 
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objects such that it vanishes for I + I' + I" = odd. For 
combinations involving an odd number (e.g. BQQ), the 
Wigner-3j symbol should be replaced with its algebraic 
approximation (B2) but with l + V + 1" = even terms set 
to zero instead. However the overall sign depends on the 
orientation of the triangle in the flat-sky approximation 
since the bispectrum is then antisymmetric to reflections 
about either axis. 



C. Weak Lensing 

In the so-called Born approximation where lensing ef- 
fects are evaluated on the the nuU-geodesics of the un- 
lensed photons, all effects can be conveniently encapsu- 
lated in the projected potential [l^jlq] 



(n) = -2 J dDg^iDM^ih),D), 



(21) 



where 



Da{D) 



-,o'£4^,.iDr (22) 



Da{D') 



For the CMB, the source distribution gs is the Thomson 
visibility and may be replaced by a delta function at the 
last scattering surface Dg = D(z ^ 10'^); for galaxy weak 
lensing this is the distance distribution of the sources. 
We explicitly relate this quantity to the more familiar 
convergence and shear in Appendix |^. Note that the 
deflection angle is given by the angular gradient = 
V</.(n).^ 

As with the temperature perturbations, we can decom- 
pose the lensing potential into multipole moments 



(n) = ^0z„,yr(A), 
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or Fourier moments as 



(n) 



(23) 



(24) 



The power spectra of the lensing potential in the all-sky 
and flat-sky cases as 



{4>%m\')) = (2^)^5(1 -l')Cjf, 



(25) 



where again C'j^^ = Cf'^. The lensing potential also de- 
velops a bispectrum in the non-linear density regime, 
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m m m ' " ' 



mm'm")) = (2vr)^-5(l - r - l")i?(T,0)" , (26) 

which is responsible for skewness in convergence maps 
and other higher order effects. Since the lensing potential 



is not affected by non-linearity until very high multipoles 
(see Fig. ^), we neglect these terms here. 

Finally, the lensing potential can also be corre- 
lated with secondary temperature and polarization 
anisotropics P,p^, so that one must also consider the 
cross power spectra 



{^lm4'l'm') — Sl^l'Sm.m'Ci 
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{X*{m\')) =(2vr)25(l- 



l')C.f* 



•'(0 



(27) 



in the all and flat sky limits. 

To calculate the power spectra of the lensing potential 
for a given cosmology one expands the gravitational po- 
tential in eqn. (|2^) in plane waves and then expanding 
the plane waves in spherical harmonics. The result is 



where 



Cr = 4^ / ^A|(fc,z)[/i™(fc)]2 



/^(fc) = j dDW'"\Dln{^^D) 
W'--{D) ^ -2F{D) D^i^^-^) 



(28) 



Da{D)Da{Ds) 



(29) 



For curved universes, replace the spherical Bessel func- 
tion with the ifltra-spherical Bessel function. In the small 
scale limit, this expression may be replaced by its equiv- 
alent Limber approximated integral jl^ ] 



(0 



27r2 



J dDDA[W'''''{D)]^Al{k,0) 
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This expression also has the useful property that its non- 
linear analogue can be calculated with the replacement 



F{DfAlik,0)~^Al{k,D), 



(30) 



where the time-dependent non-linear power spectrum is 
given by the scaling formula ||l^] and the Poisson equa- 
tion (|^). Since non-linear effects generally only appear 
at small angles, the full non-linear all-sky spectrum can 
be obtained by matching these expressions in the linear 
regime (see Fig. ||). 

Similarly, the cross correlation may be calculated for 
any secondary effect once its relation to the gravitational 
potential is known. We shall illustrate these results with 
the integrated Sachs- Wolfe effect. It contributes to tem- 
perature fluctuations as 



e^sW(n) = -2 / dD<^{^{h),D). 



(31) 



It then follows that the all-sky cross correlation is given 

by [10 

= 4-/ f A|(fc)/;-(fc)//sW(fc) , (32) 



4 



LTio-^ 




FIG. 3. Lensing power spectra. The power spectrum of the 
lensing potential is shown in the top panel as calculated by 
the fiat and all sky approaches for linear and non-linear den- 
sity perturbations. In the lower panel, the cross correlation 
with the ISW effect is shown. In both cases, a non-negligible 
fraction of the power comes from scales where the flat-sky 
approximation is inadequate. 



where 



7-ISW 



W 



[k) 

ISW(^) 



2F{D) , 



^0 



(33) 



again with the understanding that one replaces the spher- 
ical Bessel function with the ultra-spherical Bessel func- 
tions for curved universes. Similarly the flat-sky expres- 
sion becomes. 



C. 



e<j> 
(0 



k — I n ^ 



Figure || also shows the cross-correlation for the ACDM 
cosmology. 

Cross lensing-CMB bispectrum terms can also included 
but require an external measure of lensing (e.g. a galaxy 
weak lensing survey) to be observable with three-point 
correlations. 



III. FLAT-SKY POWER SPECTRA 

In this section, we calculate the effects of lensing on 
the CMB temperature (§[IIA), polarization and cross 
(§111 B) power spectra. The simplicity of the resulting 
expressions have calculational and pedagogical advan- 
tages over the traditional flat-sky correlation function 



approach [^|J^. However we also show why one cannot 
expect a flat-sky approach to be fully accurate even on 
small scales. 



A. Temperature 

Weak lensing of the CMB remaps the primary 
anisotropy according to the deflection angle 



e(n) = e(n-|- V0) 

= e(n) -I- V,0(n)V*e(n) 
1 

+ 2 



V,<?;.(n)Vj(/)(n)V'V^e(n) 



(34) 



Because surface brightness is conserved in lensing only 
changes the distribution of the anisotropics and has no 
effect on the isotropic part of the background. 

The Fourier coefficients of the lensed field then become 



9(1) 



6(1) 



dne(n)e-*'" 



e(ii)i(i,ii) 



(35) 



where 



L(l, li) = 0(1 - li) (1 - li) • li + i 1 1^0(12) (36) 

X0*(l2+ll-l) (l2-ll)(l2+ll-l)-ll. 



This determines the lensed power spectrum 



e* (1)6(1')) = (2^)2,5(1 - i')c; 



ee 



^ = (1 - PR) Cf Q + 

x[(l-ll)•ll]^ 

where 



" ^1 ^ee ^0 



(37) 



(38) 



(39) 



The second term in eqn. (^) represents a convolution 
of the power spectra. Since l^Cf"^ peaks at low Z's com- 
pared with the peaks in the CMB (see Fig. ||), it can be 
considered as a narrow window function on Cf'^ in the 
acoustic regime 200 < I < 2000. It is useful to consider 
the limit that Cf^ is slowly varying. It may then be 
evaluated at 1 — li « 1 and taken out of the integral 



(40) 



(27r)2 ' 

Note that the two terms in eqn. (|37|) cancel in this limit 
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c\ 



ee 



ee 



(41) 



This is the well known result that lensing shifts but does 
not create power on large scales. Intrinsic features with 
width A/ less than the / of the peak in I'^Cf'^ are washed 
out by the convolution (see Fig. ^). Note that in the 
ACDM model this scale is ^ ~ 40. The implication is 
that for such a model, the smoothing effect even for high 
multipoles arises from such low multipoles that the flat- 
sky approach is suspect. 

On scales small compared with the damping length 
I > 2000, there is little intrinsic power in the CMB so 
that the first term in eqn. ( ^ ) can be ignored and the 
second term behaves instead like a smoothing of Cf'^ of 
width /SI approximately the / of the peak in I'^Cf^ . Since 
Cf'^ is very smooth itself, the term is approximately, 



ee 



^ee 



2 ' 



(27r)2 



ee 



(42) 



where we have interchanged the roles of li and li — 1. The 
power generated is proportional to the lensing power at 
the same scale and may be approximated as the lensing of 
a pure temperature gradient |^] . In this limit the flat-sky 
approximation should be fully adequate. 



where recall that R was defined in eqn. (|39|). The cross 
correlations between B and O or i5 still vanish since lens- 
ing is parity conserving. Unlike the case of the temper- 
ature fiuctuations, lensing does not conserve the broad- 
band large scale power of the E and B but only the 
total polarization power. For example, lensing will cre- 
ate a B component in a field that originally had only 
an £'-component. Furthermore, lensing actually destroys 
temperature-polarization cross correlations due to the 
lack of correlation with the generated B polarization. 
From Fig. ^ one can see that the largest relative effect 
of lensing is on the correlation. 



IV. ALL-SKY POWER SPECTRA 

In this section, we treat lensing effects on the tempera- 
ture ( IV A ), polarization and cross ( IVB ) power spectra 
in a full all-sky formalism. Corrections to the flat-sky re- 
sults remain at the 10% even on small scales. Moreover, 
although the derivation appears more complicated, the 
end results for the power spectra are simple. They are 
as readily evaluated their the flat-sky counterparts and 
should be used in their stead. 



B. Polarization 



A. Temperature 



The lensing of the polarization field may be obtained 
by following the same steps as for the temperature field 

±l(n) = ±X(n + V0) (43) 
«±X(n)+V»0(n)V*±X(n) 

+ iv,0(n)V,(/)(n)V''V^±X(n), 

where we have used the shorthand notation = Q ± 
ill. The Fourier coefficients of the lensed field are then 



In the all-sky case, the Fourier harmonics are replaced 
with spherical harmonics, and the lensed field becomes 



Qim « Qhn + J dnl^"*V,</)(n)V'e(n) 

+ i j dny,"*V,0(n)Vj0(n)V^VJe(n) 



illTli /2m2 



(46) 



±X(1) = ±X(1) 



(2^)2 



,A:(li)e±2''('^'i"'^')i(l,li), (44) 



Isms 



where L was defined in eqn. (|37|). 

Recalling that ±X{1) = E{1) ± iB{l), we obtain the 
power spectra directly 

cr-{i-i'R)cr+l I ^[(i-ii)-iipC|t\i 

X [(Qf ^ + ^) + cos(4^, J(Qf ^ - ^)] , 
x[(Q^^ + Qf^)-cos(4^,J(Qf^-Q^^)], 



with the geometrical factors expressed as integrals over 
the spherical harmonics 



(47) 



jmm\m2Tn2, 
II1I2I3 



- (1 - PR) Cf"" 



J dnYi"'* (V.Fjfi) (Vj-y™^*) W^Yi:;' . 
The lensed power spectrum then becomes 



Ci^Ci+J2 ^r^^r^i + (^^"^ E , (48) 



^ = f 1 - I'R) C?"" 



xQ^^cos(2^zJ, 



[(l-li)-lipCpi^l 

(45) 



with 
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^1 = E Krr')" ' 

mim2 



(49) 



where "cc" denotes the complex conjugate and we have 
suppressed the Z-indices. 

These formidable looking expressions simplify consid- 
erably. The second term may be rewritten through inte- 
gration by parts and the identity V^Y,™ = -l{l + 



(50) 



The remaining integral may be expressed in terms of the 
Wigner-3j symbol through the general relation 



(-1) 



mi+si 



(2;i + i)(2Z2 + i)(2;3 + i) 



h h h 

Sl -S2 -S3 



where note that qY/^' 
to define 



, (51) 

—mi 1712 J 
It is therefore convenient 



hhh = ^Mh + 1) + ^3(^3 + 1) - hih + 1)] (52) 



(2;i + i)(2;2 + i)(2Z3 + i) f h h h 



47r 



Finally the Wigner-3j symbol obeys 



El h h h 
\ mi 7712 m^ 

m\m2 



h h h 
mi m2 m^ 







1 

2h + l 



Putting these relations together, we find that 



^1 = 



1 



2/4-1 



(FihhY 



(53) 



(54) 



An algebraic expression for the relevant Wigner-3j sym- 
bol is given in the Appendix. 



The second term in eqn. (48) can be simplified by re- 
expressing the gradients of the spherical harmonics with 
spin-1 spherical harmonics. As shown in Appendix the 
spin-1 harmonics are the eigenmodes of vector fields on 
the sky and naturally appear in expressions for deflection 
angles. Note that there is a general relation for raising 
and lowering the spin of a spherical harmonic []ll| , 



{l-s){l + s + l) 



s + 



1^1 



m4 



il + s){l-s + l) 



-lYr , (55) 



so that 



^yrn ^ yi([±l) [^y™^^ _ -i^rm-] . (56) 

As an aside, we note that equation (^) can alternately 
be derived from this relation and the integral (^) with 
s = ±1. 

Further, we note that spin spherical harmonics also 
obey a sum rule ||l^ 



E 



'2/ + 1 

,,Yr{n),,Yrin) = ^/ ^,,1^-^ (0) . (57) 



An 



For the spin-1 harmonics 



^iY^\0) ^ iY,-\0) 



21 + 1 
An ' 



(58) 



and the others involving si,S2 = ±1 vanish. These re- 
sults imply that 

^ V.rrV.yr* = + l)^[(m+),(m_), 

+ (m_),(m+),]. (59) 

To evaluate the second derivative term in equation (|4^), 
we again apply equation (|5^) to show that 

[(m+),(m_), + (m_),(m+),]V*V^- ,1^™ = 

-[l{l + l)-s'],Yr. (60) 



Putting these expressions together we obtain, 

S2^-hii + i)h{ii + i)^^. 



Finally combining expressions eqns. (^8|), (^4|) and (|61|), 
we have the following simple result 



(61) 



cr = [1 - i{i + m] cr + E cr^^^^^f^, m 



where 



(63) 



This expression is computationally no more involved than 
the flat-sky expression eqn. ( p8|) and has the benefit of 
being exact. Since the lensing effect even at high I in the 
CMB originates from the low order multipoles of 0, cor- 
rections due to the curvature of the sky are not confined 
to low I. We show in Fig. |l| that the correction causes 
a 10% difference in the effect. The change in Cf^ itself 
is even smaller (of order 1%). Nonetheless it is larger 
than the cosmic variance of these high multipoles and 
thus should be included in calculations for full accuracy. 
Corrections can be even larger in models with a red tilt 
n < 1 in the initial spectrum. 
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B. Polarization 

The derivation of the all-sky generalization for polar- 
ization is superficially more involved but follows the same 
steps as in the temperature case and results in expres- 
sions that are no more difficult to evaluate. The lensed 
polarization multipoles are given by 



±Xlm = ±Xlm + ^ ^ 4'limx±Xl2 



(64) 



T-rnmim2 



with the geometrical factors expressed now as integrals 
over the spin-spherical harmonics 



mmi m2m3 



r 



Noting that 



jmmim2 / -\\L j-mmirn2 



T2^nii2 



(65) 



(66) 



where L ^ I + li + I2 and recalling that ±Xim = Eim ± 
iBim, the power spectra then become 



EE 

hh 



2S1 



BB\ 



(67) 



where 



22 -S*! 
±2'S'2 



' Tmmim2\ 
.2-'lhl2 ) 



E 

mim2 

1 \ ^ jmniimim , 



(68) 



The expression for C^^ follows by interchanging EE and 
BB. The cross power spectrum is 

- cr + \ E t^r^rii + (-i)"]o25i 
iiii 

h 

with 



2 5*1 



Vlhh 2^11^12 ) 



(70) 



Just as in the case for the temperature field, these ex- 
pressions simplify considerably. The spin- 2 harmonics 
are eigenf unctions of the angular Laplacian of a tensor 



(71) 



which follows from contracting indices in equation (|& 
It then follows that 



r 

±2^llil 



= ^[^1(^1 + 1) + ^2(^2 + 1) - + 1)] 



J dn{^,Yr)Y^^ri±2K') ■ (72) 

Comparison with eqn. ( ^l| ) implies that it is convenient 
then to define the quantity 

2Fhhh = Ihih + 1) + ^3(^3 + 1) - hih + 1)] (73) 



(2Zi + l)(2Z2 + l)(2;3 + l) f h h h 

2 0-2 



so that 



2S1 



21 + 1 



47r 



i2Fihhf 



02 



Si - 



1 



21 + 1 



(-^"1(2) {2F1I1I2) ■ 



(74) 



The third term in equation ( |6q ) can be simplified by 
following the same steps for the analogous temperature 
term except for the replacement of s = with s = ±2 in 
equation (p5|). The result is 



±2'S'2 



1 2/ 

-[/(/ + i)-4];i(?i-fi)— 



2/1-1-1 



47r 



(75) 



Putting these relations together, we obtain the result for 
the power spectra 

4)i?] cr + \ E ctt^-^^^ff 



BB 



<P<P ( 2FU1I2) 



2 21 + 1 

[Qr+Qr-(-i)'(Qr'-Qr)] . 

Cf^ + 2)R] Cf ^ -t- J2 



(76) 



{Fll^l2 2E1I1I1) r^QE 



Recall that L — I + I1 + I2 and R was defined in eqn. (|62 
These expressions are plotted for the ACDM model in 
Fig. |. 
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V. FLAT AND ALL SKY BISPECTRA 



Gov = C««C««C««5,,r/;.ii5;3;;, + 5perm. 



(81) 



In this section, we consider the lensing contributions to 
CMB bispectra through the correlation with secondary 
anisotropies. We begin by reviewing the calculations 
for the temperature bispectrum as previously treated by 
We then introduce the polarization and cross bis- 
pectra which in principle have signal-to-noise advantages 
over the temperature bispectra. We illustrate the formal- 
ism with a concrete calculation of the effect due to the 
ISW secondary anisotropy. 



A. Temperature 



Contributions to the temperature bispectra from the 



cross power spectrum Cf * discussed in C| follow im- 
mediately from the first order lensing term, i.e. eqn. (^) 
for the all-sky bispectrum [0, 



B\ 



eee 



h^.ia -^'i'.'3CTQr + Sperm., 



and eqn. (B5h for the flat sky bispectrum 



B 



eee 

(Il,l2,l3 



5perm. 



(77) 



(78) 



One can show that these relations satisfy the general ex- 
pression for the correspondence between flat and all sky 
bispectra eqn. (|20|) by noting that 



I2 • h 



1 



{ll + ll - ll) 



(79) 



since the angles of a triangle is fully defined by the length 
of its sides. 

Note that there can be strong cancellation between the 
terms in the permutation in both cases. As we have seen, 
the spectrum of (j) is generally peaked to low multipoles 
implying a corresponding weighting of Cf"^ to low mul- 
tipoles for secondary anisotropies that correlate strongly 
with (j). In this case the triangles (^i,/2,^3) that con- 
tribute most strongly are highly flattened such that two 
sides nearly coincide in length h ~ h- In this case, 
contributions l\ and /| in eqn. ( |79| ) are cancelled off the 
permutation I3 <-> li leaving only a term of order t^- 

These considerations also signal problems for the flat- 
sky expressions. It is important to know what on scales 
most of the detectable signal is coming from. In the all- 
sky formalism, the signals from the m modes are added 
together with weights given by the Wigner-3j symbol 



neee 



E 



h h h 

mi 1712 TO3 



(80) 



For the small effects due to the correlation of secondary 
anisotropies with lensing, the covariance of the bispec- 
trum estimators is dominated by the Gaussian noise from 
the power spectrum H 



where the permutations are in the indices of the I' triplet. 
The overall signal-to-noise becomes 



E E«[Cov-]i3f||. (82) 

/l ^2^3 ^i^2^3 



The covariance is in general diagonal in the 6x6 blocks of 
permutations of (Zi, I2, 13) and for this simple case of the 
temperature bispectrum, the blocks are proportional to 
the trivial matrix of all ones. The result is one can take a 
simple sum over all distinct triplets or equivalently divide 
the full sum by a factor of 6, 



^ ee^ee^ee ' 



(83) 



for a cosmic variance limited experiment. For a realistic 
experiment with noise from the detectors and residual 
foregrounds, one simply replaces 



XX' 



XX' 



XX'(noisc) 



(84) 



here and below. Note that one can also construct the 
Fisher information matrix of the bispectrum along these 
lines |l9). 

Correspondingly, in the flat-sky one constructs the op- 
timal inverse- variance weighted statistic [|| (see also Ap- 
pendix ^ 



sky 



(27r) 



; ee^ee^ee ' 



(85) 



where /sky is the fraction of the sky covered. We show 
that these expressions are equivalent in the high I, /sky = 
1 limit in Appendix Thus the extra factor of /sky 
can be included in the all-sky expression to approximate 
the effects incomplete sky coverage due to exclusion of 
regions contaminated by galactic foregrounds. 

The weighting of the modes is such that the quantity of 
interest in the lensing-temperature correlation is Pcf"^ 
where the extra factor of I over the straight bispectrum 
contribution comes from the square root of the volume 
factor in Z-space. This quantity is plotted in Fig. |^ for the 
cross correlation with the ISW effect. The implication 
is that for this effect, full accuracy requires an all-sky 
approach and we shall hereafter use this to evaluate the 
signal-to-noise. 

The overall signal-to-noise as a function of the largest 
I included in the sum is shown in Fig. ^ for a cosmic 
variance limited experiment and the Planck satellite (see 
p9[ for the specification of the noise). Note that the 
Planck satellite is effectively cosmic variance limited to 
I ~ 1000 and even so the S/N is only of order a few |4|. 
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(a) Cosmic Variance 




(b) Planck 



O00 




0.1 



FIG. 4. Cumulative signal-to-noise in the bispectra as a 
function of maximum I for a cosmic variance limited ex- 
periment and for the Planck satellite. Note that for the 
cosmic variance limited case (a), bispectra involving the 
B-polarization have a substantial signal-to-noise advantage 
over the other bispectra. For the Planck satellite (b), we 
assume that the additional variance comes only from detec- 
tor noise. In practice, residual foreground contamination and 
sky-cuts to avoid them will lower the signal-to-noise further. 



B. Polarization and Cross Correlation 

Bispectra involving the E and B parity polarization 
will also receive contributions from the correlation in- 
duced by lensing. Although these signals are smaller than 
the temperature bispectrum in an absolute sense, we have 
seen that the main obstacle in detecting the temperature 
bispectrum is cosmic variance from the Gaussian contri- 
butions. 

We begin by analyzing terms that do not involve the 
B-parity polarization. For these all-sky bispectra, only 
terms with L = li + I2 + Is—even are non- vanishing, and 
we will implicitly assume that onlyeven terms are con- 
sidered. With the help of eqns. ( ^4|) and (|65|), we can 
immediately write the all and flat sky results as 



r)E0e 



Qe<t>QeE 



E<t>^e0 



-(12-13) cos 2^3lC,^*Qf^ 



(86) 



(ii • i3)(c,f + cr^^cr/) + (I2 - 13) 



where 



(87) 



The general correspondence between the flat and all 
sky expressions in eqn. ( pO| ) is established by the use of 
eqn. (|7^) the approximation discussed in the Appendix 



h h h 
2 0-2 



cos 2(^3 



h h h 




(88) 



for L =even. Th e cancellation for flattened triangles 
discussed in § V A still applies and is easiest to see in 
the flat-sky limit: the flatness of the triangles implies 

COs2(/33i ~ 1. 

For the S/N calculation, note that the covariance is 
given by 

Cov = Qf^QrQr^M'//.<i^/3.^ 



^-iBE 

-fC^^Qf ^5,^,, 5,,,/,3,, + il', <-> I',) , (89) 

so that a full calculation requires inverting this matrix 
for each distinct triplet. Since we are interested mainly 
in the order of magnitude of S/N, we can set the lower 
bound as 



> 



(T^Eee\2 

hhh '1 '2 '3 



(90) 



which amounts to ignoring duplicate triplets and replac- 
ing the remaining triplet with the average S/N of the 
set. This limit is plotted for the ISW effect in Fig. ^ as 
a function of the maximal li included in the sum. As 
expected, it is comparable to the signal-to-noise in the 
temperature bispectrum. Of course, it is experimentally 
more difficult to achieve the cosmic variance limit in the 
polarization with a realistic experiment containing detec- 
tor and foreground noise. 

There is also a qualitatively new effect from the 
polarization- lensing correlation Cf'^. However since sec- 
ondary polarization only arises from Thomson scatter- 
ing effects, we expect this contribution to be small in 
ACDM models where the optical depth during reioniza- 
tion is T < 0.3 

The EOQ bispectrum term is 



EEQ 



^E<f>^eE 

(l2-l3)(cOs2(p3lCf^q 

+ cos2^2iQrC;*^'^ 



5 



EEQ 

(Il,l2,l3 



Erjj^eE 
h 



-(ll.l2)<^C«^ + (ll-l2). 



(91) 



with covariance 



^Qf ^<5,,,. <5,,,/,3,. + {l[ - I',) , (92) 
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with which we can bound the S/N 



> 



E 



(93) 



Again, the ISW example is shown in Fig. 
Finally the EEE bispectrum is given by 



2^h^.;3C^r^r + Sperm. 



T^EEE 

-°(Il,l2,l3 



-^Ecji^EE 
h '3 

-(I2 • I3) cos2(p3i<^Qf ^ + Sperm. (94) 
with covariance 

Gov = ^Qf^Qf ^<5z,r/z.ii'5z3i. + Sperm. , (9S) 
and signal-to-noise 



c \ ^ ( uEEE\2 



N 



hhh '1 '2 '3 



(96) 



This bispectrum signal vanishes for the ISW effect. 

Bispectra involving the i3-parity polarization have dis- 
tinct properties. For terms involving one B-parity polar- 
ization term, only ^1-1-/2 + ?3=odd contributes to the 
all-sky spectrum and we implicitly assume below that 
even terms vanish. 

For the BQQ bispectrum, 

= » ( 2^w. J c^r^r + - ^3) , 

Bfi^Zis) = • I3) sin2^3iC,^'^C«^ - ih - I3) . (97) 

Again the correspondence between the flat and all sky 
expressions in eqn. ( pO| ) is established by the approxima- 
tion discussed in the Appendix 



h h h 
2 0-2 



±i sin 29331 



h h h 




(98) 



for L =odd. The sign ambiguity comes from the fact 
that a reflection of the triangle (li, I2, 13) across one of the 
axes corresponds to remappings ip n — ip or ip —ip 
and hence a reversal in sign of the flat-sky bispectrum in 
equation ( p7| ) . In this case the cancellation for flattened 
triangles discussed in § V A does not apply. However since 
sin2(/53i w 2iy93i <C 1, a suppression still exists. 
The covariance of the BQQ bispectrum is 



Gov = 5,3,- + {I', 

leading to a signal-to-noise 



E 



2(jBB(j0e(j00 

hhh '1 '2 (3 



^3) , (99) 



(100) 



In a cosmic variance limited experiment (see Fig. the 
BQQ bispectrum has signal-to-noise advantages over its 
temperature and E polarization counterparts due to the 



fact that for scalar perturbations C^^ is dominated by 
the lensing contributions themselves. Moreover, even if 
the tensor contributions are near their current limits of 
T/S ^ 0.3, the signal-to-noise is not much affected for 
I ^ 100 due to the strong damping of gravity wave contri- 
butions under the horizon scale at last scattering. How- 
ever for the Planck experiment, the detection is severely 
limited by detector noise and may also suffer further 
degradation from incomplete foreground subtraction . 
Next, the BEQ bispectrum is given by 



TjBE0 _ ■( -p ^-Bip^WK , -p ^EEf^TJcp^ 



^E4>i^0E 



EE^0<p\ 



B, 



BE0 
(li ,12,13) 



E(j)^eE 



+ sin2(^2iQf^C,'^^V 



(101) 



with a covariance 

Gov==Qf^Qf^Qf^<5,,r/i./^<5,3z. 

leading to a signal-to-noise 



C \ ^ f UB00\2 



N 



2(jBB(jEE(j00 ■ 

;i/2i3 '1 '2 ^3 



(102) 



(103) 



The signal-to-noise of this term can be greater than that 
of BQQ due to the fact that the temperature and E- 
polarization are only partially correlated in the unlensed 
sky. 
Finally, 

BF^FI = * ( 2^W2^3) <^C,f ^ + {h <-> h) , (104) 
<^b.i3) = -(I2 • l3)sin2^3i<^Qf^ - (I2 - I3) , 
with a covariance 

Gov - GF,^CF''ci''5i,v^5i,v^5i,v^ + (Z^ ^ l'^) , (lOS) 
leading to a signal-to-noise 



( r>BEE\2 



s 



N ) ^ 2C7,^^C,^^Cf^ 

'1 '2 '3 



(106) 



This signal vanishes for the ISW effect. 

Terms involving more than one B term have no contri- 
butions to first order in the correlation power spectrum. 



VI. DISCUSSION 

We have shown that a harmonic approach to weak lens- 
ing in the GMB provides a simple and exact means of cal- 
culating its effects on the temperature and polarization 
power spectra, given the power spectrum of the lensing 
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potential or convergence, and on the analogous bispec- 
tra given their power spectrum of the cross correlation 
with secondary anisotropics. Corrections to the flat-sky 
approximations appear even at high multipoles because 
even there, lensing effects arises from the large-scale fluc- 
tuations in the deflection angles. These corrections cor- 
respond to a change in the predictions at the /iK level. 
While this is a negligible change given observations to- 
day, it is above the cosmic- variance limit and should be 
included when interpreting the high-precision results ex- 
pected from Planck. 

Unlike the temperature bispectrum, bispectra involv- 
ing both the temperature and polarization multipoles of 
the CMB have the potential of producing a high signal- 
to-noise 10) detection of secondary anisotropics such 
as the ISW effects even with relatively modest angular 
resolutions I < 1000. Other secondary anisotropics such 
as the Sunyacv-Zel'dovich effect are expected to con- 
tribute even stronger signals, although their exact am- 
plitude is far more uncertain presently [Q. 

Achieving a cosmic-variance limited detection of the 
magnetic-parity polarization is a daunting challenge. 
Even signal-to-noise near unity requires detectors which 
are a factor of 3 more sensitive to polarization than those 
planned for the Planck satellite. Also of concern are the 
residual foreground contamination remaining in the maps 
after multifrequency subtraction. Our current best mod- 
els of the foregrounds indicate that with the Planck chan- 
nels and sensitivities, foregrounds and detector noise may 
enter into the polarization maps with comparable ampli- 
tudes [ pO[ . Thus improving the actual sensitivity to the 
cosmic signal beyond the specifications of the Planck ex- 
periment will not only require better detectors but also 
a better understanding of the foregrounds, perhaps with 
increased frequency coverage and sampling. 

Nonetheless, the polarization of the CMB offers the 
potential to open a new window on physical processes at 
low redshifts and the opportunity to learn more from the 
CMB than can be achieved with the next generation of 
CMB satellites. 
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APPENDIX A: ALL-SKY WEAK LENSING 
OBSERVABLES 

All weak lensing observables may be defined in terms 
of the projected potential (f> 



(n) = -2 J dDg^{D)^{^{n),D), 



(Al) 



or equivalently its multipole moments (jum in the all- 
sky formalism or Fourier coefficients (/'(I). Recall from 



eqn. ( p^ ) that is the lensing efficiency function. 

The deflection angle that a photon suffers while travel- 
ing from the source at Dg is given by the angular gradient 
of the potential a(n) = V(/)(n). Applying equation js^ ) 
to the the spherical harmonic expansion, we obtain 



Im 



E\/^^^'-[i^rm+-_iyrm-] . (A2) 



This implies that the quantity ai±ia2 is a spin ±1 object 
[ai±ia2]{n) ee ^(c ± ig),™ ^^^^"(n) 

Im 

= ±Y, VW+^)<l>irn ±iYr{h) , (A3) 



which states that the curl term cim vanishes and the gra- 
dient term 



gim = -iVKl + ^)(t>im ■ (A4) 
The power spectrum of the angular deflection is then 

{9l'm'9lm) = ^l.l'^m.m'Cf^ 

= Sl,l,Sm,m'l{l + l)Cf'^, (A5) 

with the curl power vanishing. This accounts for the 
factors of l{l + 1) in equations involving the angular de- 
flection [e.g. eqn. (|6^)]. 

The corresponding flat-sky quantity is given by the 
decomposition [see eqn. (C8)] 



[ai ± iQ!2](n) = ±i 



dH 



[c±zg](l)e=^*('^'-'^)e''-", (A6) 



with c(l) ~ and 



5(1) = -im) 



(A7) 



These relations also give the bispectrum of the deflec- 
tion angle in terms of bispectrum of the lensing potential 
in the obvious manner. 

The convergence (k) and shear (71,72) are familiar 
weak lensing observables from galaxy weak lensing stud- 
ies [|l^ . Although they are not directly needed for CMB 
studies, they are of interest for cross-correlation of galaxy 
weak-lensing maps and the CMB. An equivalent all-sky 
lensing treatment is given by [|2^ . 

These quantities are given by the second derivatives 

ViVjCj) = ngij + (71 + 472) (m+ (g) m+),y 

+ (71 - «72)(m_ ® m_)y , (A8) 

convergence where gij is the metric on the sphere. 

For the all-sky harmonics, it is useful to note that equa- 
tion (pa) implies 
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and hence 



1(1 + 1) 



Y, 



1 {1 + 2] 



^Y[^{vn+ ® m+) + _2>"r(m- ® m 



(A9) 



<n)^~J2^lil + l)^i^Yrin), 



Im 



7i 



(n) ±i72(n) = XI 9 



1 G + 2)! 



f-2V (Z-2)! 

Consequently, the power spectra are related as 
P{l + iy 



±^Yr in). (AlO) 



cr 



Cf 



4 

1 (/ + 2)! 

4(^-2)!"'' ' 

-h{l + l)C^\ 



2 V (/ - 2) 



(All) 



where the e shear power spectra is defined in the same 
way as that of the E polarization and X ~ ld,E, B. The 
/? shear power is the analogue of the B polarization power 
and vanishes for weak lensing. 

In the fiat-sky limit, these expressions become 



n{n) = - 
71 (n) ± 172 (n) = - 
so that 



1 f <F\ 



2 J (27r)2 
1 [ d^l 



f(l>{\)e 



(A12) 



{2ny 



±2i{ipi-ip) ^il-n 



^(0 



(1) 



(A13) 



These relations also give the bispectrum of the shear and 
convergence in terms of the bispectrum of the lensing 
potential 

Khis = iMh + mih + i)hih + i)KZ, , 



_ 1 l {h + 2)\{h + 2)l{h + 2)l ^^^ 
^hhh - 8 Y (/i - 2)! ih - 2)1 ik - 2)! '^'^'3 ' ^A^*^ 

with a similar relation for the flat-sky bispectra. 

APPENDIX B: WIGNER-3J EVALUATION 



1. Exact Expressions 

The expressions for the power spectrum of the lensed 
temperature and polarization distributions involve spe- 
cific sets of Wigner-3j symbols that can be efficiently 
evaluated. The expression for the temperature involves, 
a set which has a closed algebraic form: 



h h h 




{L-2h)\{L-2h)\{L-2h)\^^'^ 



(Bl) 



(i + 1)! 

for even L = li + I2 + and zero for odd L. 

The required set for the polarization does not have an 
exact closed form expression. However it may be equally 
efficiently evaluated for our purposes with the realization 
that in the sums, we require 



^1 h h 
mi 1712 



(B2) 



for fixed /2, ^3, wi, TO2, TO3 and all allowed li. The recur- 
sion relations for the Wigner-3j symbol, 

WAi^+iwi^+i + Bi^wi^ + (h + l)Ai^wi^-i = , (B3) 

where 



A, = ^il - [h - h)^^{i2 + h + 1)2 - il^q - ml , 

Bl, = -{2h + l)Ml2 + l)mi - kik + l)mi 

-;i(/i + l)(m3-m2)], (B4) 

allow us to generate the whole set at once . For a sta- 
ble recursion, one begins at the minimum and maximum 
1 1 values 



^imin = max(|;2 - h]), \mi\) , 

^Imax = ^2 + ^3 , 



(B5) 



with ifij^jj,^ = wi-,_^^^ = 1 and carries the recursion in both 
directions to the midpoint Zimid in the range (or any non- 
vanishing entry in the vicinity). One then renormalizes 
either the left or right recursion to make the wi-,^^-^^ agree. 
The remaining overall normalization is fixed by requiring 



Y^{2h + l)wl^l. 
h 



and 



sgn(w,i„^J = (-1) 



h-ls—mi 



(B6) 



(B7) 



Putting these relations together, we obtain the full set of 
symbols as required. 
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FIG. 5. Wigner-3j function and approximation. An ex- 
ample of the Wigner-3j symbol relevant to the polarization 
calculation with I2 = 100, m2 — 0, I3 = 50, ma — —2 is shown 
as calculated from the recursion relations (solid upper) and 
analytic approximation (dashed). The difference is shown be- 
low (solid lower). 



2. Approximations 

We can use the general relation between the all and 
flat sky bispectra of equation (pO|) compared with the 

to 



explicit calculation of the flat sky bispectrum in § V B 



develop an high-Z approximation for the speciflc symbol 
in the polarization calculations. The comparison implies 
that 



for L 



Then 



h h h 
2 0-2 

I1 + I2 + h 

cos 2<y931 = 



cos2v.3i(^o' 0' 0' j , (B8) 
even. By the law of cosines, 

(B9) 



2;2 

3 



III 



2 



\L/2 



2 Z-^/g 



(L/2)! 



{Ll2^h)\{Ll2-l2)\{Ll2~h)\ 



{L-2h)\{L-2h)\{L-2h)\ 



1/2 



for L = even. For odd values of L, we use the relation 



W h h 
2 0-2 



±isin2(y93i 



h h h 




(BIO) 



and fix the overall sign ambiguity by an explicit evalua- 
tion. By the triangle relations. 



sin 2(^31 = T^mL ~ 2h){L ~ 2l2){L ^ 2h)]^/^ 

z| — — z| 



Putting this together with equation (B2) and fixing the 
sign ambiguity, we obtain 



h h h 
2 0-2 



(-1) 



L-l 1 fll- If 

2 



2 ;2 
3 



72 72 

'l''3 



(Bll) 



(i/2)! 



{L/2-hy.iL/2-i2y.iL/2-hy. 

L{L-2li){L-2l2){L~2k) 
(i-2/i)!(L-2/2)!(i-2Z3)!li/2 



(i + 1)! 

for L = odd. The half integer factorials are defined by the 
gamma function x\ = T(l + x). By explicit calculation 
we find that these expressions are valid to better than 3% 
of the rms amplitude of the symbol when averaged over 
neighboring / for all li — Ih — hl ^ 25 and ^2 + ^3 — ^ 25, 
i.e. for triangles that are sufficiently far from being flat. 
Near zero crossings, the fractional error can be large but 
the absolute error remains a small fraction of the rms. A 
typical case is shown in Fig. ^. 

These relations may be useful in cases where only a sin- 
gle symbol is needed. However for the lensing calculation 
where the whole set is required, the recursion relations 
are as efficient as the approximation and are exact. 



APPENDIX C: FLAT AND ALL SKY 
CORRESPONDENCE 

1. Harmonics 

We establish here the correspondence between the all 
and flat sky harmonic cocflicients of spin zero (scalar), 
spin one (vector) and spin two (tensor) quantities on the 
sky. 



Following 1 23 1 , let us begin by introducing the following 
weighted sum over the multipole moments of the field 
X = E, B, or (f> for a given I and its inverse relation. 



Xil) 



An 



21 



IS- 



(CI) 



The goal is then to show that this quantity is the Fourier 
coefficient of the flat-sky expansion. 

Spin-0 quantities, such as the temperature flucutations 
and the lensing potential, are decomposed as 



X{h)=Y^Xi^Yr{n). 



(C2) 
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For small angles around the pole, the spherical harmonics 
may be approximated asj^] 



(C3) 



V ZTT 

and the expansion of the plane wave 

gii-n ^ ^i'"j^(;6i)e™(v-v,) 



I 

Thus 

Im 



I e , 



(C4) 



(C5) 



which is the desired correspondence. 

Spin-1 quantities like the deflection angles are decom- 
posed as 



(C6) 



Here one notes that 



^,Yr^±Je'^"'{^.±^^y)Yr, (C7) 



and thus 



Irn 



(2^ 



-X{\)e 



±i(ipi-ip) a-: 



(C8) 



Finally, spin-2 quantities like the polarization are de- 
composed as 



(C9) 



Ira 



Here one notes that 



1 



±2Yr^i^e^^'^{d.±^dyfYr, (CIO) 



*Note that our definition of F;'" differs from the usual one 
by ( — 1)™ to conform with the spin spherical harmonic con- 
vention [G2l. 



and thus 



Irn 



2 il-n 



dH 

(27r)2 



,X(l)e±2''('^'-'^)e*'-", (Cll) 



as desired. 



2. Power Spectra 

The correspondence between power spectra then fol- 
lows from the relationship between the harmonics 



xJ(l-l'). 



,imipi 



dipi>e 



(C12) 

We then expand the delta function in plane waves func- 
tions 

dn 



5(1-1') 



(2^)2 



27r f dn 



(C13) 



■m—rii \rm ^\rm irmpi—im ^p,/ 
I lit lid 



Integrating over the azimuthal angles Lpi , (fv collapses the 
sum to 



XX' 



C^f J dnYf"^*Y,7"^' 

XX' 



which proves the desired relation in cqn. (po|), 

r<xx' ^ r^xx' 



(C14) 



(C15) 



3. Bispectra 

The correspondence between bispectra is established in 
exactly the same way as with the power spectra. The only 
difference is that the expansion of 5(1' — 1) in eqn. (C14) 
is replaced with that of d{li 4- 12 + I3) leading to 



I I' I" 



m m m' 



I I ^ii'i 



xx'x'- 



(C16) 



BfJ',^; I dnY-^Yi"^ Yi,"^ 



= B 



xx'x" 



I V l"\ f I V I 



{i,i',i") I 



(2^-H)(2r-H)(2/" + l) 
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This establishes the relation 



TjXX'X'' 

^ii'i" 



I V I" 


xx'x" 



[21 + l){2l' + l)(2l" + I) 



in 



B 



(1,1', I") 



(C17) 



Note that we have implicitly assumed that the bispec- 
trum only depends on the the magnitudes (Z, Z', I") so that 
it may be removed from the azimuthal integrals. This is 
not true for terms not involving the magnetic parity. In 
this case, the sign of the flat-sky bispectrum depends 
on orientation but we find empirically that a similar re- 
lationship holds up to a sign ambiguity as discussed in 



for the flat sky bispectrum The extra factor of (27r)^ 
compared with the power spectrum is from the extra 
delta function in the noise term. One can show that 
these expressions agree in the high-Z limit by restoring 
the integration over Z3, expanding the delta function into 



spherical harmonics as in eqn. (C14), and integrating over 
azimuthal angles, 



h 



lidli 



idln 



+ I3) 
'(27r)5 



hdli / Z2CZZ2 / hdl 



h h h 




(C23) 



4. Signal-to-Noise 

Here we establish the correspondence between the all 
and flat sky signal-to-noise statistics for the case of di- 
agonal contributions to the covariance matrix [Gov = 
diag(Var)], 



^ Var 

Im 



= ^(2Z + l 



Var 



(C18) 



For the flat sky case, one defines a weighted sum of 
Fourier harmonics 

P = J d^lW{l)X{l)X{-\), (C19) 

with optimal weights given by W{1) = C^-^^ /Va.i- from 
which one calculates the signal-to-noise (P)^ / {P'^) as 



Var 



2/s 



ky 



Idl^^^^. 
Var 



(C20) 



where we have used the fact that (5(0) « V/{2tt'^) = 
/sky/TT. These expressions agree in the high /-limit and 
imply the familiar result that /sky should multiply the 
signal-to-noise of angular power spectrum measurements 
given incomplete sky coverage. 

The bispectrum signal-to-noise similarly is 



I1I2I3 



^XXX\2 

Var ^ 



XXX\2 



h h h 
mi 1712 TO3 



Var 



(C21) 



for the all sky bispectrum and 



jXXX 



With the general correspondence of bispectra from 
eqn. (C17), this becomes 



S 



n) -^f^^yj dhj dhj dk 



(nXXX\2 



(C24) 



which proves the equivalence of the signal-to-noise for 
high-/ and /^ky = 1. 
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